Fermi's golden rule applied to the 7 decay in the quasicontinuum of 46 Ti 
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Particle-y coincidences from the 46 Ti(p,//y) 46 Ti inelastic scattering reaction with 15-MeV protons are uti- 
lized to obtain y-ray spectra as a function of excitation energy. The rich data set allows analysis of the coinci- 
dence data with various gates on excitation energy. For many independent data sets, this enables a simultaneous 
extraction of level density and radiative strength function (RSF). The results are consistent with one common 
level density. The data seem to exhibit a universal RSF as the deduced RSFs from different excitation energies 
show only small fluctuations provided that only excitation energies above 3 MeV are taken into account. If 
transitions to well-separated low-energy levels are included, the deduced RSF may change by a factor of 2 — 3, 
which might be expected because the involved Porter-Thomas fluctuations. 

PACS numbers: 21.10.Ma, 25.20.Lj, 27.40.+Z, 25.40.Hs 



I. INTRODUCTION 

Fermi's golden rule predicts the transition rate from one 
state to a set of final states in a quantum system. The the- 
oretical foundation, which has been successfully applied in 
many disciplines of physics, was first established by Dirac in 
1927 0] and emphasized by Fermi in his book in 1950 0. 
The rule is based on first-order perturbation theory, where 
the transition matrix element is assumed to be small. In nu- 
clear physics, this assumption is well fulfilled for j5 and y de- 
cay. Thus, we take the validity of the Fermi's golden rule as 
granted, rather than testing this rule. In this work, we study 
the y decay between states in the quasi-continuum of 46 Ti and 
apply Fermi's golden rule to disentangle the y strength and 
level density. 

The Oslo nuclear physics group has developed a method to 
determine simultaneously the level density and the radiative 
strength function (RSF) from particle-y coincidences. These 
quantities provide information on the average properties of ex- 
cited nuclei and are indispensible in nuclear reaction theories 
as they are the only quantities needed for complete description 
of the y decay at higher excitation energies. 

The nuclear level density can be determined reliably up to a 
few MeV of excitation energy from the counting of low-lying 
discrete known levels [3|. There is also reliable level-density 
information from neutron resonances at higher excitation en- 
ergies; however, this information is restricted in energy as well 
as spin range. 

The most rich experimental information on the RSF was ob- 
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tained from the study of photonuclear cross-sections JU, and 
thus limited only to energies above the particle separation en- 
ergy. It was established that the giant electric dipole resonance 
(GEDR) dominates the RSF in all nuclei. The information on 
RSF below particle separation energy is significantly less. It 
has been obtained mainly from the Oslo method and (n,y) and 
(Y>Y) reactions. 

The Oslo method, which is applicable for excitation ener- 
gies below the particle separation, is described in detail in 
Ref. J3). In this work, we report on results obtained for 
the 46 Ti nucleus, which has two protons and four neutrons 
outside the doubly-magic 40 Ca core. The advantages of the 
(p,p') reaction compared to the commonly used ( 3 He, 3 He') 
and ( 3 He, 4 He) reactions are much higher cross-sections and 
better particle resolution. This allows us to make a detailed 
study of y-decay as a function of excitation energy. 

It can be discussed if the concept of one unique RSF is 
valid for light nuclei and in particular at low excitation ener- 
gies. This question together with the applicability of the Oslo 
method for light systems as titanium is the main subject of this 
work. 

In Sec. II the experimental results are described. The nu- 
clear level density and RSF are extracted in Sec. Ill, and in 
Sec. IV, the application of Fermi's golden rule and the Brink 
hypothesis is discussed. Summary and conclusions are given 
in Sec. V. 



n. EXPERIMENTAL RESULTS 

The experiment was conducted at the Oslo Cyclotron Labo- 
ratory (OCL) with a 15-MeV proton beam bombarding a self- 
supporting target of 46 Ti with thickness of 1.8 mg/cm 2 . The 
target was enriched to 86% 46 Ti with 48 Ti (11%) as the main 
impurity. This small admixture of 48 Ti is not expected to play 
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FIG. 1 : (Color online) Singles (a) and coincidence (b) proton spectra 
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FIG. 2: (Color online) The particle- y coincidence matrix for 46 Ti. 
The y-ray spectra have been unfolded with the Nal response func- 
tions. 



a major role. This is supported by the fact that at low exci- 
tation energies where the level density is small, we could not 
identify any important contributions of 48 Ti into the y-spectra. 
In addition, both titaniums are expected to behave similarly. 

Particle-y coincidences were measured with eight Si AE — 
E particle telescopes and the CACTUS multidetector sys- 
tem 0. The Si detectors were placed in forward direction, 
45° relative to the beam axis. The front (AE) and end (E) de- 
tectors had a thickness of 140 jim and 1500 jJ.m, respectively. 
The CACTUS array consists of 28 collimated 5" x 5" Nal(Tl) 
y detectors with a total efficiency of 15.2% at Ey = 1 .33 MeV. 
The singles -proton spectrum and protons in coincidence with 
y-rays are shown in Fig.[T] 

In total, 110 million coincidence events were collected in 
one week with a beam current of 1 .5 nA. Using reaction kine- 
matics, the measured proton energy was transformed into ex- 
citation energy of the residual nucleus. In this way, a set of 
y-ray spectra is assigned to a specific initial excitation energy 
Ej in 46 Ti. Furthermore, the y-ray spectra are corrected for 
the known response functions of the CACTUS array following 
the procedure described in Ref. [7 |. The unfolded coincidence 
matrix (E^Ey) of 46 Ti is shown in Fig. [2] 

The coincidence matrix displays vertical lines that represent 
yrast transitions from the last steps in the y-cascades. How- 
ever, there are also clear diagonal lines where the y energy 
matches the direct decay to the ground state (£, = Ey) and to 
the first and second excited states at 889 keV (2 + ) and 2010 
keV (4 + ), respectively. These y-rays are primary transitions 
in the y-cascades. By studying the energy distribution of all 
primary y-rays originating from various excitation energies, 



information on the level density and RSF can be extracted si- 
multaneously. 

An iterative subtraction technique has been developed to 
separate out the first-generation (primary) y-transitions from 
the total cascade [8|. The subtraction technique is based on 
the assumption that the decay pattern is the same whether the 
levels were initiated directly by the nuclear reaction or by y 
decay from higher-lying states. This assumption is automat- 
ically fulfilled when states have the same relative probability 
to be populated by the two processes, since y-branching ra- 
tios are properties of the levels themselves. If the excitation 
bins contain many levels, it is likely to find the same y-energy 
distribution from this set of levels independent of the type of 
population. However, the assumption is more problematic if 
the decay involves only a few (but not one) levels within the 
energy bin. 

Fermi's golden rule predicts that the decay probability may 
be factorized into the transition matrix element between the 
initial and final states, and the density at the final state 0: 

^ f = T \(f\H'\i)\ 2 p f . (1) 

Realizing that the first generation y-ray spectra P(Ei,Ey) 
are proportional to the decay probability from £, — > Et = 
Ei — Ey, i.e. Xi-jff, we may write the equivalent expression of 
Eq. ([T]l as: 

P(Ei,Ey) oc ,^ fPf , (2) 

where is the y-ray transmission coefficient, and P/ = 

p(Ei — Ey) is the level density at the excitation energy Ef af- 
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ter the primary y-ray emission. This expression does not allow 
us to extract simultaneously ,%^f and p from the experimen- 
tal P(Ej,Ey) matrix. To do that, either one of the factorial 
functions must be known, or some restrictions have to be in- 
troduced. According to the Brink hypothesis [9|, the y-ray 
transmission coefficient is independent of excitation energy; 
only the transitional energy Ey plays a role. Thus, we replace 
&U f with 3{E y ), giving 

P{E u E y ) «* ST(E y )p h (3) 

which permits a simultaneous extraction of the two multi- 
plicative functions. We then fit about N 1 j2 data points of the 
P matrix with 2N free parameters. A least % 2 fit is then pos- 
sible, since it is many more data points than fit parameters; in 
the present case we have 150 free parameters to fit 2240 data 
points. 

At low excitation energy, the y decay is highly dependent 
on the individual initial and final state; therefore, we have ex- 
cluded y-ray spectra originating from excitation energy bins 
below £, =5.5 Me V. 

It is well known that the Brink hypothesis is violated when 
high temperatures and/or spins are involved, see Ref. iflOl and 
references therein. However, in the present Oslo experiment, 
the temperature reached is low (T ~ 1.5 MeV) and is assumed 
to be rather constant. The dependency on spin is of minor im- 
portance. The measured ratios of the y feeding into the ground 
band indicate a low spin window of / ~ — 6Ti. At excita- 
tion energy E — 8 MeV the ratios are 57:100:9 for the 2+, 4+ 
and 6+ states, and at E ~ 10 MeV they are 55:100:10, which 
are equal within the error bars for extracting these ratios. Of 
course, at low excitation energy the spin distribution fluctuates 
due to a few levels (and spins) present within each 118 keV 
excitation bin. 

In principle, the Brink hypothesis assumed in expression <|3j 
could bias the analysis, so that the Oslo method in the next 
turn validates the Brink hypothesis. This issue will be ad- 
dressed in Sec. IV, where we find one common level density 
(according to Fermi's golden rule), which in turn results in 
one universal RSF in the quasi-continuum of 46 Ti. 



III. LEVEL DENSITY AND RADIATIVE STRENGTH 
FUNCTION 

In our first investigations, we rely on the Brink hypothe- 
sis from the expression Q and factorize the first-generation 
y-matrix P(Ej,Ey) into one transmission coefficient 2?(Ey) 
and one level density p(E). Since the decay between states at 
low excitation energy cannot be treated within a statistical ap- 
proach, a cut of the matrix with Ey> 1.8 MeV and 5.5 MeV 
< Ei < 10.0 MeV was used to exclude clear non-statistical 
decay routes 1 . The least % 2 fitting of the two multiplicative 
functions follows the iterative procedure of Ref. 151 . 



1 The lower excitation cut concerns only the initial excitation energy £, ; one 
still needs y-spectra originating from excitation regions down to the ground 
state in order to subtract higher-order generations of y-rays. 



To demonstrate how well the procedure works, we compare 
in Fig.|3]for some initial excitation energies £, the experimen- 
tal first-generation spectra P with the ones obtained by multi- 
plying the extracted 3T and p functions. The agreement be- 
tween calculated and experimental first-generation spectra is 
excellent for decay from higher excitations; however there are 
locally strong deviations where the calculated spectra fall out- 
side of the experimental error bars for populations of lower ex- 
citations energies, as seen in the £, = 5.6 and 6.4 MeV gates. 
These deviations could be the consequence of Porter- Thomas 
fluctuations [11], which will be further discussed in Sec. IV. 
The general good agreement holds also for all the other 40 
spectra (not shown) included in the global fit with the same 
S?(Ey) and p(E) functions. 

The experimental statistical errors are very small as seen 
in Fig. [3] Thus, the deviations are due to other sources of 
errors. For example for Ej = 6.4 MeV, the 3" ■ p prediction 
overestimates around Ey = 3 MeV and underestimates around 
Ey = 5.5 MeV with several standard deviations. Thus, there 
are indications that a common 3 and p function that simulta- 
neously fits the P(Ej,Ey) matrix, could not be found. The sys- 
tematic errors behind these deviations could be due to several 
factors as experimental shortcomings, assumptions behind the 
Oslo method, the Brink hypothesis |9| and, most probable, 
Porter- Thomas fluctuations. Keeping all these possibilities in 
mind, the results of Fig.|3]are very gratifying. 

There exist infinitely many 3" and p functions making 
identical fits to the data as the examples shown in Fig. [3] 
These functions can be generated by the transformations 

p(Ei-Ey) = Aexp[a(£,--£ r )]p(£,-£ y ), (4) 
&{E y ) = Bexp(aEy)3 r {E r ). (5) 

In the following, we will try to determine the parameters A, a, 
and B. This information is not available from our experiment 
and has to be determined from other experimental results. 

First, the level density at high excitation energy is nor- 
malized to the level density at the neutron separation energy 
p(S„). This data point is calculated from neutron resonance 
spacings Do (see, e.g., Ref. [ 12 1) with a spin distribution given 
by Gl 
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2/+1 
2a 2 



exp[-(/+l/2) 2 /2(J 2 ], 



(6) 



where E is excitation energy and / is spin. 

There exists no neutron resonance data for 46 Ti. We there- 
fore estimate p(S„) from the parameterizations of von Egidy 
and Bucurescu |[T4l using the back-shifted Fermi gas (BSFG) 
model, which reads 



Pbsfg(£) = rj 



exp(2-\ 



12A'W 4 ff' 



(7) 



where a is the level density parameter, U = E — El is the in- 
trinsic excitation energy, and El is the back-shift energy pa- 
rameter. The spin cut-off parameter a is given by lfl4l 



a 2 = 0.0146A 5/3 



1 + ^1+4-aU 
2a 



(8) 
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FIG. 3: (Color online) Comparison of experimental first-generation spectra (squares) and the ones obtained from multiplying the extracted S7 
and p functions (red line). The initial excitation energy bins E, are 118 keV broad. The error bars represent the experimental statistical errors. 



A being the nuclear mass number. 

Figure |4] shows the extracted total level densities for tita- 
nium isotopes with known resonance spacings Do at S„ (filled 
triangles). The resonance spacings for I = neutrons only 
give the level densities for one or two spins and only one par- 
ity. In order to extract the level density for all spin and parities, 
we use Eq. (|6]i and assume equally many positive and negative 
parity states at S„. This is supported by combinatorial quasi- 
particle models ITsHlTl . where the numbers of positive and 
negative parity states at S n are predicted to be the same. 

The points connected with lines are based on the semi- 
empirical estimate of von Egidy and Bucurescu [14| with a 
common scaling of 77 = 0.5 in order to match qualitatively 
the experimental p (S„) points. The estimated value for 46 Ti is 
p = (4650 ± 1000) MeV- 1 at S„ = 13.189 MeV (see Fig.Q. 
The error bar chosen reflects roughly the general deviation 
between the global estimates and the points derived from neu- 
tron resonance data. 

Now the scaling (A) and slope (a) parameters of the level 
density can be determined as shown in Fig. [5] The normal- 
ization of the level density is fitted to the known level density 
around 3.5 MeV of excitation energy, and to the extrapola- 
tion from p (S n ) using the BSFG model with parameters sum- 
marized in Table [i] By choosing other p(S„) values (within 
the assumed uncertainty of ±1000 MeV -1 ), the logarithmi- 
cal slope of the p and & will change accordingly, as the a- 
parameter of Eqs. Q [5]) has to be adjusted. The p curve is 
well fixed at ~ 3.5 MeV and will rotate around this point with 
different choices of a. Level densities in the 14 — 19 MeV 
excitation region of 46 Ti have been measured from Ericson 
fluctuations, see Ref. |[T8l and references therein. However, 
the values reported by the various experimental groups differ 
within a factor of ten. The measurement of the Ohio group 



lfT8l is shown in Fig. [5] for comparison (see triangle point at 
15.5 MeV). 

One should stress that the level density found from the ex- 
periment is based on the spin and parity levels populated in the 
(p,p') reaction. Thus, the normalization to the total level den- 
sity described above, rests on the assumption that the struc- 
ture of the level density remains approximately the same if all 
spins and parities are included. This assumption is reasonable 
fulfilled according to Ref. ifTBI . where the level density for 
spin windows of 2 — 6h and — 30h have been calculated for 
46 Ti within a combinatorial quasiparticle model 1171 . From 
these estimates, our measurement includes 70 — 80% of the 
total level density and the level density fine structures for the 
two spin windows are very similar. 

The measured level density describes nicely the known 
level density up to around 4 MeV of excitation energy. The ex- 
perimental resolution of p is about 0.3 MeV at low excitation 
energy, as seen for the 2 + state at 889 keV. At around 3 MeV, 
we see an abrupt increase in level density due to the breaking 
of proton and/or neutron pairs coupled in time-reversed or- 
bitals (Cooper pairs |fl9l ). This results in ~ 10 times more lev- 
els, making it difficult to determine the level density at higher 
excitation energies with conventional spectroscopic methods. 

It remains to determine the scaling parameter B of the trans- 
mission coefficient &(Ey), Here we use the radiative width 
(Ty) at S„ assuming that the y-decay is dominated by dipole 
transitions. For initial spin / and parity n, the width is given 

by ma 

xp(S n -E r ,I f ), (9) 
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FIG. 4: Deduced total p(S„) from neutron resonance spacings (trian- 
gles). The data point for 46 Ti (square) is estimated by extrapolations 
from the BSFG model with global parameterization of von Egidy and 
Bucurescu 1141 (cirles with lines to guide the eye). 



where the summation and integration run over all final levels 
with spin If that are accessible by y radiation with energy Ey 
and multipolarity El or Ml. Further details on the normaliza- 
tion procedure are found in Refs. ll5l l2TI . 

Since there exists no neutron resonance data, we again have 
to rely on systematics from other isotopes. From Fig. |6| we 
estimate an average 7- width of (Ty) = (1200 ±500) meV at 
S n . The large uncertainty in (r^) gives an absolute normal- 
ization of S?(Ey) and the RSF within about ±50%. However, 
this uncertainty does not influence the energy dependence. 

The deduced RSF for dipole radiation can be calculated 
from the normalized transmission coefficient ST(JLy) by El 



f(Ey) = 



1 gXgy) 
2k El 



(10) 



The normalized RSF is shown in Fig. [7] For comparison, the 
GEDR data from Ref. Il23ll are also shown, which have been 
translated from photo neutron cross section a to RSF by 
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3n 2 h 2 c 2 



(11) 



J 7 



Unfortunately, there is a large energy gap between our data 
ending at a E y = 10 MeV and the GEDR data that start at 
14 MeV. 

Our data on the RSF display a minimum near 4 — 6 MeV 
and some structures at lower y-ray energies. The uncertain- 



FIG. 5: Normalization of the nuclear level density (filled squares) of 
46 Ti. At low excitation energies, the data are normalized (between 
the arrows) to known discrete levels (solid line). At higher excitation 
energies, the data are normalized to the BSFG level density (dashed 
line) going through the point p(S„) (open square), which is estimated 
from the systematics of Fig. [4] For comparison a data point 1181 
obtained from Ericson fluctuations are shown (black triangle). 



ties in the value of p at S„ will change the slope of the RSF, 
and thus also the degree of low-energy enhancement. How- 
ever, the lines of Fig. [7] show that the enhancement is not very 
sensitive to reasonable choices of the value p(S„). 

Such or similar enhancement has been seen in several 
light nuclei with mass A < 100, see Ref. [24 1 and references 
therein. There is still no theoretical explanation for this very 
interesting phenomenon. 



TABLE I: Parameters used for the extraction of level density and 
radiative strength function in 46 Ti. 



S n a E\ 0" 
(MeV) (MeV" 1 ) (MeV) 


P (s„) <r r (s„)) 

(MeV" 1 ) (meV) 


13.189 4.7 -1.0 4.0 


4650(1000) 1200(500) 



IV. FERMI'S GOLDEN RULE AND BRINK HYPOTHESIS 

According to Fermi's golden rule, it is possible to factor- 
ize out the level density from the y-decay probability. Since 
the level density is a unique property of the nucleus, the same 
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extracted level density is expected from the decay probability 
deduced at different excitation energy regions within the ex- 
perimental restrictions, in particular the energy resolution of 
CACTUS. In order to see if the Oslo method gives a unique 
level density, we have divided the data set into three statis- 
tically independent initial excitation-energy regions, namely 
E t = 5.5 - 7.0 MeV, 7.0 - 8.5 MeV and 8.5 - 10.0 MeV, and 
applied the same methodology as described in Sec. III. The 
results shown in Fig. [8] are very satisfactory. The different 
data sets give approximately the same level density. Thus, the 
disentanglement of level density and transmission coefficient 
from y-particle coincidences seems to work very well accord- 
ing to Eq. Q. 

The fact that we measure approximately the same level den- 
sity function for primary y-ray spectra taken at different ex- 
citation regions indicates that the RSF depends only on y- 
ray energy and not on excitation energy. This seems to in- 
dicate the validity of Brink hypothesis, which will be tested 
more thoroughly in the following. In principle, the test could 
be performed by dividing the data set into even more initial 
excitation-energy regions. However, the statistics of the ex- 
periment does not permit such an approach, and a different 
approach has been used. 

By accepting the level density obtained with the global fit 
of all relevant data (see lowest curve in Fig. [8]), we may further 
investigate the transmission coefficient in detail, and thus the 
validity of the Brink hypothesis. We first adopt the solutions 
3~ and p from Sec. Ill and write 

JS(E,)P{Ei,E r ) « &{Ey)p{Ei-Ey). (12) 

The normalization factor for each initial excitation bin is de- 
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FIG. 7: (Color online) Experimental radiative strength function for 
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(3650 MeV -1 ) level density at S n are also shown. For comparison 
the GEDR data from the (y.abs) reaction |23 | are shown (triangles). 



fined by 



J^d£ y ^(£ y )p(£,--£ y 
J^dE r P{Ei,E r ) 



(13) 



The degree of correctness of the approximation ( 12 1 is typi- 
cally illustrated by the fits in Fig. [3] 

In the following, we will investigate the dependency of the 
deduced RSF on initial and final excitation energy. If the 
Brink hypothesis is valid, there exists the same RSF for all ex- 
citation energies in this nucleus. We will call it the universal 
RSF. In reality the Porter- Thomas fluctuations involved influ- 
ence the RSF obtained from different excitation regions. We 
will use the term deduced RSF in the following for the quan- 
tity obtained from experimental data. The deduced RSFs are 
expected to fluctuate around the universal RSF and the fluc- 
tuations are expected to be stronger if the number of transi- 
tions used in the determination of the deduced RSF is smaller. 
Thus, the deduced RSFs extracted from data sets involving 
initial and final regions of high level density should be much 
closer to the universal RSF than the RSFs deduced from re- 
gions of low level density. 

Since there exists only one unique level density, we con- 
struct the counterpart to Eq. ( 12 1 in the case that the transmis- 



sion coefficient depends on the initial excitation energy: 

^K'(Ei)P(Ei,E r ) a & (Ey, E()p (Ej —E r ), (14) 



7 



10< 



10 s 



Q. 

>» 
*^ 

gio 2 

T3 

> 
0) 



10 









1 


5.5-7.0 MeV 






7.0-5 


1.5 MeV 




X27 


8 T 5-10.0 MeV 






j$> 5.5-10.0 MeV 














- , ■ , 1 


1 < 1 1 1 1 < 1 1 





2 4 6 8 

Excitation energy E (MeV) 



10 
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upper curves). The lower curve is the result for the whole energy 
region and is identical to the one of Fig. [5] 



where Jf' is determined analogously to Eq. ( 13 1. We expect 
that ,^{Ei,Ey) fluctuates on the average around 3T(Ey). Thus, 
it is reasonable to expect that jY' ~ jY , which gives a trans- 
mission coefficient of: 



P{Ei,E y ) 
p(E,-E y 



(15) 



Similarly, the transmission coefficient as a function of the final 
excitation energy Ef = Ej — Ey is given by 



&(Ef,Ej)mjr(Ef + E, 



■Ey,Ey) 



PiEf) 



(16) 



The validity of the approximations ( fl4} and ( p"5j ) is demon- 
strated in Fig. [9] by comparing the deduced RSF from this 
approximation (lines) with the independent fits of and p 
(data points). The RSFs determined for the whole energy re- 
gion 5.5 — 10 MeV are very similar, especially for low Ey, 
as shown in the lower panels. Also the similarity in the de- 
tailed structures of the two methods is recognized, although 
some differences are present. The deviations are largest for 
the high y-energy part of the RSFs populating the lowest + , 
2 + , and 4 + states, where large Porter-Thomas fluctuations 
are expected (these fluctuations are not included in the error 
bars). However, the overall good resemblance encourages us 
to study the detailed evolution of the RSFs as a function of 
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FIG. 9: (Color online) Gamma-ray strength functions extracted from 
various initial excitation bins E,. The data points of the three 
upper deduced RSFs are from statistically independent data sets. 
The data points in the lower panel are identical with the RSF of 
Fig. [7] The RSFs displayed as red lines, are evaluated from the ra- 
tio P(Ej,Ey) I p(Ef) (see text). The resemblance between the two 
methods confirms that the level density is common for the various 
excitation regions, in accordance with Fermi's golden rule. 



( 15 1. Since the approximation ( 16 1 is a simple transformation 
using Ef = Ej—Ey, we may also investigate the dependencies 
of the RSF on final excitation energy Ef. 

In Figs. [TO] and [IT] the RSFs f(Ei,E Y ) and f(E f ,E Y ) are 
shown for eight excitation regions. For all these deduced 
RSFs, we use one common level density in the evaluation 



of the approximations ( 15 1 and (16 1. The data points of the 
experimental P matrix cover only a certain region in E, and 
E Y , making restrictions on the deduced RSFs. Thus, f(Ej,Ey) 
is limited to 5.5 MeV < Ej < S„ and 1.8 MeV < E y < Ej. 
The limits for f(E f ,E r ) are < E f < S„ - 1.8 MeV and 
1.8 MeV < Ey < S„ — Ef. For a consistency check, we have 
tested that the average RSF for all Ef energies equals the one 
for all Ej energies (not shown here). 



initial excitation energies using the approximations ( 14 1 and 



Figure 10 shows that the deduced RSF changes as a func- 
tion of E(, when low excitation energy is populated after the 7- 
emission. The lower, left panel, where two consecutive RSFs 
are plotted together, illustrates this. Here, the bumps seen at 
high y-ray energies are the artifacts of the decay to specific 
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FIG. 10: (Color online) Deduced RSFs from various initial excitation bins E,. The RSFs are evaluated from the ratio P(Ej,Ey)/ p(Ef) as 
described in the text. 



isolated levels below 3 MeV of excitation. This is also the 
reason why apparently the fluctuations are typically a factor 
of 2 — 3 in the panel where all deduced RSFs are plotted to- 
gether (lower, right panel). In the plot of all deduced RSFs, 
we see a minimum at about 4 — 6 MeV and some interest- 
ing structures at low y-ray energies. These structures and the 
minimum are independent of the uncertainty in the log-slope 
of the level density, as shown in Fig. [7] 

The various f(Ef,Ey) plots in Fig.[TT|are difficult to com- 
pare due to different limits appearing at both low and high 
y-energies for the various Et regions. For example, the first 
three spectra do not reveal the region of low y-energy en- 
hancement because Ey > 5.5 MeV — £> > 3.7 MeV. These 



/ 



spectra represent the decay to the + , 2 + and 4 + ground band 
states, respectively, where the experimental energy resolution 
makes some overlap between these states. In general the var- 
ious spectra show strong fluctuations in the deduced RSFs 
when the y emission ends up at low excitation energy, typi- 
cally Ef < 3 MeV. The panel with all deduced RSFs plotted 
together (lower, right panel) shows approximately the same 



scattering of data points as in Fig. 10 



It is thus clear that the experimentally deduced RSFs for 
which states below 3 MeV of excitation energy are involved 
in the y decay are very different. Figure [5] shows that this 
excitation region coincides with a region of few and well sep- 
arated low-lying levels with specific structures. The bumps in 
the deduced RSFs are specific to the low energy level scheme, 



and the changes between the various deduced RSFs are large 
due to Porter- Thomas fluctuations. 

In order to show the similarity of the deduced RSFs in the 
case of strongly suppressed Porter-Thomas fluctuations, we 
have compared the two uppermost excitation energy gates E, 
and Ef in the lower, left panels of Figs. [lOland 11 respec- 



tively. The deduced RSFs are here extracted for ydecay be- 
tween states in quasi-continuum, except for the data points 
with Ey > 7 MeV in Fig. 10 where the final excitation en- 



ergy is < 3 MeV. The deduced RSFs extracted from the quasi- 
continuum region behave similar to the deduced RSF from the 
all-over fit displayed in Fig. [7] 

The good agreement between the deduced RSFs at higher 
energies is consistent with the expectation of suppressed 
Porter- Thomas fluctuations due to more initial and final lev- 
els in the evaluation of the y strength. These results strongly 
indicate that the concept of an RSF, which is independent on 
excitation energy, is valid already at relatively low excitation 
energies in 46 Ti. At the same time, the results give us con- 
fidence that Oslo method works reasonably. The differences 
in deduced RSFs from lower energies indicate that the Porter- 
Thomas fluctuations are so poorly suppressed that it is difficult 
to predict the shape of the universal RSF. This is not very sur- 
prising and the differences in the deduced RSFs seem to be 
consistent with the expected fluctuations. 

In order to display the small differences in the deduced 
RSFs obtained from regions of higher level density, we have 
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FIG. 11: (Color online) Deduced RSFs populating various final excitation bins Ef. The RSFs are evaluated from the ratio P(Ef + 
Ey,Ey)/p(Ef) as described in the text. It should be noted that there are no final states at Ef = 1.2— 1.8 MeV, however the experimental 
resolution (see Fig.pl is responsible for including the 2 + and 4 + ground band states in this gate. 



taken the four highest gates shown in Fig. 10 and only con- 
sidered 7 energies up to 5.1 MeV. Thus, these statistically in- 
dependent RSFs are evaluated in quasi-continuum with initial 
and final excitation energies of roughly E, = 8 — 10 MeV and 
Ef = 3 — 7 MeV, respectively. The deduced RSFs are pre- 



sented in the upper panel of Fig. 12 



For this data set we evaluated the relative deviations by 



f.j-ifi) 
ifi) 



(17) 



where the index i represents the y energy and j is the initial 
excitation energy. The average strength at each y energy is 
estimated from the four individual RFSs 



ifi) 



(18) 



In the lower panel of Fig. 12 the r i; values are plotted showing 



the relative fluctuations from the mean value at each 7 energy. 
The average ratio for all data points («; = 28 and n,- = 4) is 
taken as 



—IK 

if 



(19) 



giving r ~ 6%. The differences in the deduced RSFs can eas- 
ily be interpreted as remnants of the Porter-Thomas fluctua- 
tions. However, a quantitative estimate of the fluctuations is 



difficult to determine since the experimental statistical errors 
(see error bars) are also of the same order. 



V. CONCLUSIONS 

The level density and radiative strength function for 46 Ti 
have been determined using the Oslo method. Similar level 
density functions have been extracted from statistically inde- 
pendent data sets covering different excitation energies. This 
gives confidence to the Oslo method, since the disentangle- 
ment of the level density by Fermi's golden rule predicts one 
and only one unique level density, independent of the data set. 

The deduced RSF displays an enhancement at low y-ray en- 
ergy where we see a bump around 3 MeV and another struc- 
ture at energies near 2 MeV. A similar enhancement has been 
seen in several other light mass nuclei and is still not ac- 
counted for by present theories. 

A method to study the evolution of the deduced RSFs as a 
function of initial and final energy regions has been described. 
The deduced RSFs are found to display strong variations for 
different initial and final excitation energies if transitions to 
the lowest excitations are involved. The reason for the violent 
fluctuations of a factor of 2 — 3 is that only a few isolated lev- 
els are present at low excitation energies with Ef < 3 MeV. 
The differences in the RSFs obtained from a few number of 
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transitions, that can be explained as a consequence of Porter- 
Thomas fluctuations of individual intensities, show that this 
energy region cannot be used for determination of the univer- 
sal RSF. Even though, the deduced RSFs based on a restricted 
number of transitions still indicate that the decay is governed 
by a universal RSF. 



However, the present work shows that it is possible to get 
more precise experimental information on the universal RSF. 
By imposing restrictions on the initial and final excitation 
energies, the RSFs for the decay between states in quasi- 
continuum can be extracted (i.e. for Ef > 3 MeV). The results 
from this selected data set show that the decay is consistent 
with an RSF which is idependent of excitation energy within 
less than ~ 6 % already at these relatively low excitations. 



In summary, the observation of almost identical level den- 
sities and RSFs extracted from different data sets in quasi- 
continuum, indicates that the Oslo method works well. Pro- 
vided that we use data from the quasi-continuum, a universal 
RSF in the light mass region of the 46 Ti nucleus can be ex- 
tracted. 
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